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Abstract 



This paper considers the mean variance portfolio management problem. We ex- 
amine portfolios which contain both primary and derivative securities. The challenge 
Qh ■ in this context is due to portfolio's nonlinearities. The delta-gamma approximation is 

£H ' employed to overcome it. Thus, the optimization problem is reduced to a well posed 

quadratic program. The methodology developed in this paper can be also applied to 
pricing and hedging in incomplete markets. 
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\D \ 1 Introduction 

<^> . The main objective in portfolio management is the tradeoff between risk and return. 

Markovitz, [8] and [9] studied the problem of maximizing portfolio expected return for a 
given level of risk, or equivalently minimize risk for a given amount of expected return. 
One limitation of Markovitz's model is that it considers portfolios of primary assets only. 
^ \ Recent works looked at the optimal management of portfolios containing primary and 

derivative assets. Here we mention [10] and [2]. In [10| . the author introduces a technique 
for optimizing CVaR (conditional value at risk) of a portfolio. The paper [2] notices 
that the problem of minimizing CVaR for a portfolio of derivative securities is ill-posed. 
Furthermore, [5] shows that this predicament can be overcome by including transaction 
costs. 

There are some papers which consider portfolio optimization with non-standard asset 
classes; we recall [T], [3], and [7]. In a continuous time model pQ looks at the problem of 
maximizing expected exponential utility of terminal wealth, by trading a static position in 
derivative securities and a dynamic position in stocks. In a one period model [I] analyses 
the optimal investment and equilibrium pricing of primary and derivative instruments. 
The paper [7J shows how to approximate a dynamic position in options by a static one 
and this is done by minimizing the mean-squared error. 
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valuable advice, suggestions and a thorough reading of the first version. 
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To the best of our knowledge this paper is the first work to consider the mean variance 
Markovitz portfolio management problem in one period model with derivative assets. For 
a portfolio containing many assets (primary and derivatives) the estimation of the correla- 
tion matrix is a big challenge. Practitioners solve this difficulty by projecting portfolios on 
a smaller numbers of factors. If parametric approaches are used (we work in a multivariate 
normally distributed returns framework), the projection method creates another problem, 
since the projections are often nonlinear; in order to overcome it the delta gamma approx- 
imation is employed. The delta-gamma approximation is well-known and often used in 
risk management and portfolio hedging. In the industry practice this approximation does 
well for small time intervals. By performing the delta gamma approximation the portfolio 
management problem is reduced to a quadratic program. Another challenge may come 
from covariance matrix of the factors not being positive definite. This issue appeared in 
some financial optimization problems; e.g., for portfolios of stocks the sample correlation 
matrix is just an approximate correlation (because is usually built from inconsistent data 
sets) and hence not positive definite. This problem is addressed by [3] and [6]. These 
works focus on the extraction of a positive semi-definite variance-covariance matrix, ob- 
tained through the solution of a second-order conic mathematical programming problem. 
It is a way to convexify an a priori non convex problem. In [3] and [6], the smallest distor- 
tion of the original matrix which satisfies the desired properties (e.g. being a correlation 
matrix) is obtained by using Frobenious norm. 

The results of our paper can be applied to the problem of pricing and hedging in 
incomplete markets. For instance we can consider instruments written on nontradable 
factors (e.g. temperature) and they can be hedged with tradable instruments which are 
highly correlated (this procedure is called cross hedging). Take as an example weather 
derivatives (e.g. HDD or CDD); energy prices are considered as the traded correlated 
instrument (in California a high correlation can be observed between temperature and 
energy prices). Perfect hedging is not possible in this paradigm. Minimizing the variance 
of the hedging error can be captured as a special case of mean variance optimization 
problem for a portfolio of primary and derivative instruments. A survey paper on mean- 
variance hedging and mean- variance portfolio selection is 

Another possible application of our results is the hedging of long maturity instruments 
with short maturities ones. As it is well known, the market for long maturity instruments 
is illiquid, thus the issuers use (static) hedging portfolios of the more liquid short maturity 
instruments. The interested reader can find out more about this in [S]. 

The paper is organized as follows: Section 2 presents the model. Section 3 introduces 
the delta gamma approximation. Section 4 presents the reduction to quadratic programs. 
Section 5 is an application to pricing and hedging in incomplete markets. 



2 The Model 

Portfolios returns are derived from the return of individual positions. In practice, it is 
not good to model the positions individually because of their correlations. If we have m 
instruments in our portfolio we would need m separate volatilities plus data on Tn ("~ 1 ) 
correlations, so in total m ^ +1 ) pieces of information. This is hard to get for large m. 

The resolution is to map our m instruments onto a smaller number of n risk factors. 
The mapping can be nonlinear (e.g. BS (Black Scholes formula) for option). Let us assume 
that the factors are represented by a stochastic vector process S = (Si, S2, • ■ • ,S n ), which 
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at all times t £ (0, oo) is assumed to be of the form 



S t = fit + ZW t . (2.1) 

Here /i is the vector of returns, S is the variance-covariance matrix which is assumed 
positive definite, and Wt is a standard Brownian motion on a canonical probability space 
(tt,J-t,J-). The value of portfolio at time t, denoted V(S,t), is of the form 

m 

V(S,t) = J2 x k(t)V k (S,t), (2.2) 

k=l 

where Vk(S, t), k = 1, • • • m, represents the value of the individual instruments (mapped 
onto the risk factors), and Xk(t), k = 1, • • • m stands for the number of shares of instrument 
k held in the portfolio at time t. We choose the portfolio mix Xk(t), k = 1, • • • m such that 
the portfolio return AT/ over time interval [t, t + At] 

AV = V{S + As,t + At) -V(S,t), (2.3) 

is optimized in a way which is described below. It turns out to be more convenient to work 
with the vector of actual proportions of wealth invested in the different assets. Thus, at 
time t € (0, oo), we introduce the portfolio weights Wk{t), k = 1, • • • m, by 

In the following we posit the Markowitz mean-variance type problem; given some exoge- 
nous benchmark return r e (t), at time t an investor wants to choose among all portfolios 
having the same return r e (t), the one that has the minimal variance Var(A^) : 



(PI) min Var(AF) 

w 

such that E(AV) = r e (t), 

m 

k=l 

Another possible portfolio management problem is to choose the portfolio with the minimal 
variance: 

(P2) min Var(AF) 

w 

m 
k=l 

There are some difficulties in solving (PI) and (P2). First, we might be short of AV^ 
moments information. Because AV nonlinearly depends on the change of factors, it is not 
obvious what distribution AV would follow even if we perfectly learn the p.d.f of AS. The 
situation would not get much better if we only require the moment information of AV. 
The integration for moments might be still hard to calculate explicitly. One way out of 
this predicament is to use delta gamma approximation. 
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3 Delta-Gamma Approximation 



The delta-gamma approximation states that a portfolio change during a given time period 
resulting from the change of underlying factors can be approximated by some second order 
polynomial function, the coefficients of which are given by the portfolio's sensitivities such 
as the delta, gamma and theta. It is an important tool in risk management and hedging; 
for instance, to hedge a portfolio of derivatives with respect to the underlying's change, 
the delta-gamma approximation is employed to match sensitivities of the portfolio with 
those of the hedge instruments. 

Mathematically speaking, this approximation is a second order Taylor expansion of the 
portfolio change AV over the time interval [t, t + At] : 

AV fa 5V = ^-At + 5 T AS + IaS t TAS, (3.1) 
at 2 



where 



Since 



then 



dV d 2 V . _ 



V(S,t)=J2 x k(t)V k (S,t), 

k=l 



^ = f^ = X>*(*)<t ^ := fl ' i = V"-n, k = l,---m, (3.2) 
* fc=i % 

d 2 V m d 2 V 

T v = dsM = ^ Xk{t)T ^ r '/ : ^ 

1 3 k=l % 3 

It is well known that this approximation performs well as long as the time interval At is 
not too big. At this point we formulate the approximated versions of (PI) and (P2), as 
follows: 

(P3) min Var((5F) 

such that E(6V) = r e (t), 



k=l 



(P4) min V&i(5V) 

w 

m 

Y,w k (t)V k (S,t) = l. 



k=l 



The next step is to reduce (P3) and (P4) to quadratic programs and this is done in the 
next section. 
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4 Quadratic Programs 

Let us first consider the case of one asset, m = 1. In the light of (I2.ip . AS" ~ AA(/x, HyJ At). 
For computational convenience we assume \x is the zero vector and At = 1. Next, replace 
the vector of correlated normals, AS, with the vector of independent normals Z ~ A/"(0, /). 
This is done by setting 

AS = CZ with CC T = S. 
In terms of Z, the quadratic approximation of Ay becomes 

AV & 5V = a + (C T 5) T Z + ^Z T {C T TC)Z, 

with 

dV 

a = -At (4.1) 

At this point it is convenient to choose the matrix C to diagonalize the quadratic term in 
the above expression and this is done as follows. Let C be a square matrix such that 

CC T = S (4.2) 

(e.g., the one given by the Cholesky factorization). The matrix ^C T TC is symmetric and 
thus admits the representation 

X -C T YC = UAU T , (4.3) 

where A = diag(Ai, • • • , A n ), and U is an orthogonal matrix such that UU T = I. Next, set 
C = CU and observe that 

CC T = CUU T C T = S, (4.4) 
hj T YC = hj T (C T TC)U = U T {UAU T )U = A. 



Thus, with 
we get 



b = C T 5, (4.5) 



AV ss 5V = a + b T Z + Z T AZ := Y. 



4.1 Moment Generating Function 

In this subsection, we explore the moment generating function of Y and further derive the 
mean and variance of Y. In the light of 

n 

Y = J2( X i Z i+ b i Z i)+ a ( 4 - 6 ) 

i=l 

n / h \ 2 n hi 



2XJ ^ i\ 

1 = 1 v ' 1=1 



it follows that the random variable Y is student distributed, being (up to a constant) the 
sum of squared independent normally distributed random variables. Thus, it is well known 
that 

E(0Y) = expfa(0)), (4.8) 
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where 



7 ? (0) = afl + ^(0) = a0 + ^- \ l-29X- ~ l0g(1 
j= i j=i V i 



29X 



•3) > 



(4.9) 



3=1 3- 

for all 9 satisfying max,- 9Xj < i. Direct computations lead to 



d(e"W) 



r/0 



exp(? ? (6>))^ 



exp(7?(6 1 )) 



+ £ 



;V l-20A J )-0 2 (-2A i ; 
(1 - 2^Aj) 2 



-2Ao 



1 - 20A,- 



and 

d 2 (e^>) 
d9 2 



exp(r](9)) 



lJ^/29b 2 (l-29X j )-9 2 (-2\ j ] 

fl ^L — 



(1 - 29Xjf 



1 - 29 Xj 



-2\j 



l 2Aj -(1-20A,) 2 " 



(26 2 - 806 2 Aj + 40&|Aj)(l - 2#Aj) 2 - (29b 2 {l - 29Xj) + 2 6 2 2A j )((-2A i )2(l - 26 Xj)) 

(1 - 29Xj)* 

Thus, the first and second moments of Y are 



and 



Hence, 



E(y 2 



E(y) 

d 2 ( e "W) 



dB 



it 



d# 2 



= a + Aj 
=o i=i 



i=i j=i 



Var(Y) = E(y 2 ) - E 2 (Y) = Yjtfj + 2A|). 

i=i 

In order to ease the notations we assume that V(S, t) = 1, so the vector of shares x equals 
the vector of proportions w (also notice that for simplicity we dropped the t dependence 
of w). We would like to express the mean and variance of Y in terms of x. In the light of 
f|4.3l) . (14. 2p . ()3.3p and trace properties it follows that 



E(Y) = a + ^2X j 

3=1 



a + tr (UAU T ) 



(4.10) 



a + -tr [C 1 TC 

1 / m 
a+ 2 tr ( XJ x i r i S 



T 

a + £ p, 
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where the vector p is defined by 



P - 2 



^(tr^.tr (r 2 S),...,tr(r m S)) T . 
As for the variance, recall that with b of (j4.5|) it follows that (see (j3.2[) and (j4.4|) ) 

n 

J2 b l = b T b = {C T 5) T C T 5 = 5 T C T C5 = -x T tx, (4.11) 

k=l 

where 

t = 2Af T SM and M = (M^) = (<5*), i = l,...,n, k = 1, . . . , m. 

The matrix £ is positive semidefinite. In the light of (|4.4|) and trace properties it follows 
that 

fe=i 4 

= -tr (rcc T rcc T ) 

= ^tr(rsrs) 



4 
1 

4 

where the matrix Q is defined by 



HjfH) 



i =1 i^i 

37 X tj 



Qij = tr (PEFE) , i = 1, ■ ■ ■ m, j = !,-■■ m. (4.12) 
Therefore, we end up with 

m 1 

Var(y) = *£(b 2 k + 2Xt) = ^ T (t + Q)x. (4.13) 

j=k 

Thus, from (|4.1U|) and (|4.13|) . the portfolio problem (P3) (recall that x = w) becomes 

(P5) min -x T (t + Q)x 
x 2 

s.t. a + x p = r e , 



^ Vfc(i, S)x k = 1, 



k=l 



and (P4) becomes 
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(P6) min — x T (T, + Q)x 
x 2 

m 

s.t. ^2v k (t, S)x k = 1. 

k=l 

It turns out that the problem (P5) has a similar form with the classical mean variance 
portfolio problem, a quadratic objective function and linear constrains. Notice that the 
matrix S + Q is positive definite. This comes from + Q)x = Yax{5V) > 0. We wrap 

up our findings in the following Theorem. 

Theorem 4.1. (P3) is equivalent to (P5), which is a convex quadratic program, and thus 
solvable in polynomial time. 

5 Quadratic Hedging 

The results we established so far can also be applied to hedging. The motivation comes 
from incomplete markets. Indeed, financial markets are fundamentally incomplete. It is 
well known that in incomplete markets perfect hedging is not possible. One way to solve 
this problem is to consider quadratic hedging; that is, minimize the variance of the hedging 
error. Let F be a payoff of the form F = Vi(St+Atit + At), for some map V±. We would 
like to hedge this payoff by some instruments which are of the form V k (S, t), k = 2, • • • , I 
(with I possible less than n, whence the incompleteness). For simplicity assume that in 
this market borrowing and lending of cash is done at zero interest rate (this can be easily 
achieved if one takes the zero coupon bonds as numeraire). Given the number of shares 
(xi, X2, ■ ■ ■ , xi) in the hedging portfolio, the hedging error is 

l+l 
k=l 

with x\ = —1, and A VJ+i (S, t) = 1. Therefore, the problem of minimizing the variance 
of hedging error is of the form (P2). The initial amount needed to finance the hedging 
portfolio is 

X l+1 + V!(St,t). 
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